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Preface : 



This report develops an analysis of the transient response 
of a gun barrel subject to repeated firing. A finite-element 
approach is used. A piecewise cubic approximation of the de- 
flection ensures continuity of deflections and slopes at the 
net points. Initially, linear and angular displacements and 
velocities are prescribed. This arbitrariness in the initial 
conditions allows for residual motions from the preceding shot. 
A digital computer program based upon the analysis has been de- 
veloped. The program may be used for the period during which 
the projectile leaves the muzzle by setting terms pertaining to 
the projectile and the gas pressure equal to zero. If desired, 
the program may be used to determine the static deflection of 
the gun barrel. Horizontal transverse vibrations of the barrel 
may be treated by discarding the acceleration due to gravity. 
The projectile is treated as a point mass. Consequently, pro- 
jectile balloting is not included in the theory. The effects 
of a gun-barrel tuning mass, eccentricity of the breech, 
eccentricity of the recoil mechanism, stiffness and damping of 
the supports, and prescribed motion of the gun foundation, as 
well as Bourdon pressure and projectile-barrel friction are 
included in the model. In the computer program, the time in- 
tegration routine employs Newmark's beta method. The computer 
program listing is not included here. However, information on 
the computer program is available by writing to 
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Professor A. P. Boresi 

Department of Theoretical and Applied Mechanics 
University of Illinois 
Urbana, Illinois 61801 
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NOTATIONS 



t = time. Dot denotes time derivative 

x = axial distance along the barrel from the breech 
(Fig. 1) 

y = deflection of the barrel; y = 3y/3t 

0 = notation of cross section of barrel. If shear 

deformation is discarded, 0 = 3y/3x = y f 

L = length of barrel 

a = radius of the bore 

n = coordinate of the tuning mass 

m^ = mass of the tuning mass (Fig. 2) 

I T = moment of inertia of the tuning mass about a diametral 
axis through its center of mass (Fig. 2) 

m = mass of the projectile 

M = mass of the breech block 

M = mass of the barrel 

m = mass of the part of the barrel in the range greater 

than x (Eq. 26) 

T = moment of inertia of the breech block about a trans- 
verse axis through its center mass 

g = acceleration of gravity 

E = Young ' s modulus 

e,d = lengths defining the eccentricity of the breech 
(Fig. 1) 

a = angle that the barrel forms with the horizontal 
(Fig. 1) 

oig = spring constant for the recoil mechanism (assumed 
constant 

Pg = viscous damping coefficient for the recoil mechanism 

v = v(t) = speed of the projectile relative to the breech 

? = <-(t) = displacement of the projective relative to the 

breech. v = £ 
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F 



u = 
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x 



1 ’ 
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A . 
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a ^ 



a . 
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v 3 

P . 
J 
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B is 

A a6 



X T+1 

e ' , e" 
P 

a = 



p(t) = gas pressure behind the projectile 

F(t) = axial frictional force of the projectile on 
the barrel 

u(t) = recoil displacement of the barrel (positive 
backwards, Fig. 1) 

P(t) = tension in the barrel at section x due to 
inertia (Eq. 26) 

mass density of barrel 

x , , = net points. x, = 0, x = L 

2’ ’ n+1 ^ 1 ’ n+1 

x . - x.. This is the length of the j-th interval 

J + 1 J 

y(xj, t) 

constant for translational spring support at point x. 

J 

constant for rotational spring support at section x. 

constant for translation damper at point x^ 

constant for rotational damper at point x ^ 

value of P at x = xj . Note that P . depends on the 

location E, of the projectile '- 1 

0.(t) = value of 9 at x = x . 

J J 

shear coefficient (Eq. 8) 

coefficients for the Bourdon effect (Eqs. 24 and 25) 

coefficients for the axial inertia effect (Eqs. 28 
and 31) 

net point ahead of the tuning mass. x T < n < x_, , 

(Fig. 2) 1 ” 1+1 

= eccentricities of breech spring and damper (Fig. 1) 

p(x) = mass of barrel per unit length 

number of segments; n+1 = number of net points 
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1. INTRODUCTION 



A finite-element method is used to analyze the transient 

response of a gun barrel due to repeated firing. The dynamic 

behavior of the barrel may be computed for the period in which 

the projectile is in the barrel; of particular importance are 

the slope and the lateral velocity at the muzzle. By setting 

to zero terms pertaining to the projectile and the gas 

pressure, one may follow the gun response also for the period 

after the projectile leaves the muzzle. The barrel is modeled 

as a tapered elastic beam with finite degrees of freedom by 

subdividing it into n intervals (finite elements) with n+1 

net points (nodes). The generalized coordinates are the 

deflections y. and the rotation 9. of the barrel cross section 
3 3 

at the nodes. Hence, there are 2(n+l) generalized coordinates. 
Two additional coordinates that enter the equations are the 
axial recoil displacement u of the barrel and the axial dis- 
placement Z, of the projectile relative to the breech. However, 
these coordinates do not enlarge the system of differential 
equations if the axial movements of the barrel and the projec- 
tile are regarded as known functions of time. 

The breech is treated as a rigid mass with its center 
offset from the axis of the barrel. Viscous dampers (dashpots) 
may be located at all nodes or they may be provided at only a 
few points, since in the computer program, a damper can be 
eliminated by setting its coefficient equal to zero. Also the 
offset of the breech may be set equal to zero if the center 
of mass of the breech lies on the axis of the barrel. The 
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structure supporting the barrel and the breech is modeled as 
a set of independent translational and rotational springs and 
dashpots at the net points. At nodes where there is no such 
support , the spring constants and/or dashpot constants are 
set equal to zero. A tuning mass may be located at an arbitrary 
point along the barrel. 

The foundation upon which the gun barrel and its supports 
rest is regarded as a rigid body that moves in a prescribed 
way. In general, it has six degrees of freedom. However, in 
a linear theory, oscillations of the barrel in a vertical plane 
are uncoupled from transverse horizontal displacements. Also, 
the rotation of the barrel about its axis is uncoupled. Con- 
sequently, attention may be restricted to a plane motion of the 
barrel. Oscillations of the barrel in a fixed vertical plane 
are considered. Horizontal lateral oscillations may be super- 
imposed. The program for horizontal lateral oscillations is 
the same as for oscillations in a vertical plane, except that 
effects of gravity are eliminated. 

The differential equations of motion may be obtained from 
Hamilton's principle, as the Lagrangian equations for a noncon- 
servative system. In terms of the virtual work <5W and the 
kinetic energy T, Hamilton's principle states that 



where t- denotes time. However, because certain contributions 



t 



1 




(oT + 6W)dt = 0 



( 1 ) 
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to the differential equations of motion can be written down 
directly in terms of well-known finite element procedures, it 
is expedient to do so. Hence, in this report Hamilton's 
principle is used mainly as a guide to verify the form of 
certain unusual effects such as the Bourdon effect, the axial 
inertia effect and other pressure effects that are inherent in 
gun systems. 

Contributions to 6W come from several sources. One source 
is the action of gravity on the breech, the barrel, the 
projectile, and the tuning mass. Another source is the axial 
frictional force of the projectile on the barrel. It introduces 
a nonconservative effect. Another contribution comes from the 
so-called '’Bourdon effect", which arises because, in the bent 
barrel, the area of the bore above the neutral plane, on which 
gas pressure acts, is slightly greater than the area below the 
neutral plane. Still another effect contributing to 5W is the 
axial inertia of the recoiling barrel, which initially exerts 
a straightening influence. The contribution of the strain 
energy of the barrel due to flexural stresses, etc. may be 
introduced as a negative contribution to 6W. This latter effect 
essentially results in a finite element model of the gun barrel 
in the final differential equation form for the gun system. 
Accordingly, a finite element beam model is used to incorporate 
this effect in the differential equations, rather than derive 
this effect from Eq. (1) directly. Finally, to account for the 
effects of the supporting structure, the support structure is 
regarded as sets of linear rotational and translational decoupled 
springs and dashpots located at the nodes of the finite element 
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model of the barrel. These elements are assumed to be attached 
to the gun foundation which may undergo prescribed motions. 

The spring constants and the damping constants for the trans- 
lational and rotational springs and dashpots at node j are 

respectively (a., a'.) and ( y • , u^). At a node where there is 
J J 3 3 

actually no supporting spring or dashpot, the spring constants 
and/or dashpot constants are set equal to zero in the computer 
program . 

In general, for a finite element model of the barrel, the 
net virtual work may be expressed as 

n+1 n+1 

6W = Q . 6y . + ^ Q 5 0 . + R3u + S3? (2) 

Z — i J J Z — t j j 

j=l j=l 



where n = the number of elements (n+1 = number of nodes) and 

where Q., Q'. , R’, and S are the components of generalized force; 
3 3 

i.e., the coefficients of the generalized coordinates u, Z, , y., 

J 

9. in the expression for 6W. 

The kinetic energy T is the sum of the kinetic energies of 

the breech, the projectile, the barrel, and the tuning mass. 

It is a function of the generalized coordinates u, £, y., 0., 

J 3 

where dots denote derivatives with respect to time t. Conse- 
quently, the differential equations (the Lagrangian equations) 
of motion are, by Hamilton’s principle, 



_d_ 3T 3T 
dt " 3y j 
j = 1, n+1. 




_d_ 3T 
dt 3 0 



ae j = S 



(3) 
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( 3 )cont ' d 



_d_ 3T _ 3T = R _d_ 3T _ 3T 

dt 3u ' ' ’ dt 



where Qj , , R, S are the components of generalized force (Eq. 2). 

The Lagrange equations for u and £ are irrelevant, if we know 
u(t) and £ ( t ) as functions of time t. It is to be noted that 
3T/3yj = 3T/30J = 0, since T = T(u, £, y y 9 j ) . 

By the usual procedures of finite element methods, Eqs . (3) 
may be written in the form 



Mx + Cx + Kx = f(t) 



(4) 



where M = M(t), C = C(t) and K = K(t), are the mass, damping 
and stiffness matrices for the gun system, f(t) is the driving 
force vector and x is the vector of generalized coordinates, 



= (y, 



n + 1 



9 n + l } 



(5) 



The Lagrange equations (Eqs. 3 or 4) are linear second order 

differential equations for the unknown functions y.(t) and 9.(t). 

j 3 

They are adaptable to various integration methods, such as the 
Runge-Kutta method, Newmark's beta method, Houbolt's method, 
etc. Discontinuities occur at the instantaneous location of 
the projectile, since, for example, the gas pressure behind the 
projectile, causing the Bourdon effect, terminates there. 

However, discontinuities do not obstruct the finite element 
method. This is an advantage over the direct approach via the 
differential equations of beams, since with the latter formula- 
tion, the action of the projectile is described with a Dirac 
delta function or a Heaviside step function. These difficulties 
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are avoided by the finite-element treatment. A computer program 
based on Newmark's beta method has been devised for solving the 
differential equations for arbitrary initial values y.(0), 0.(0) 

J J 

and y . ( 0) , 8.(0). Vertical bending and horizontal bending of 
J J 

the gun barrel are decoupled in linear theory. Both of these 
types of deformation are covered by the program. For horizontal 
motion, there is no gravitational effect. 

The program provides the history of the motion of the barrel 
from the initial instant at which the projectile is fired, until 
the projectile leaves the muzzle. In addition, it may be used 
for the period after the projectile has left the muzzle. To 
account for this latter period, all terms pertaining to the 
projectile (e.g., gas pressure, projectile mass, friction) are 
dropped . 

To implement the program, we must know the recoil displace- 
ment u(t), the projectile displacement £(t) along the barrel, 
the base pressure p(t) on the projectile, and the frictional 
force F(t), between the barrel and the projectile. If £(t) and 
p(t) are known, F(t) is determined by Newton's law applied to 
the axial motion of the projectile. In applying the program 
to an actual gun system, the spring constants (a^ , a \ ) and the 

damping constants (y., y '. ) of the supporting structure must be 

J J 

estimated. Because of the complicated nature of the structure 

and the peculiar damping devices in use, this may be difficult. 

Of particular significance for the accuracy of shooting are 

the muzzle displacement, y , , , muzzle rotation 0 ., and the 

n+l n+1 

muzzle velocity y n+ ^ at the instant that the projectile leaves 
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the barrel. The program permits a study of these quantities 
corresponding to various initial motions and deflections of 
the barrel at the instant of firing. 

2. DESCRIPTION OF THE MATHEMATICAL MODEL 

Figure 1 represents a gun barrel with mass M and elevation 
angle a. Axes (x, y) move axially with the barrel as it recoils, 
but they undergo no lateral displacement. The projectile is 
regarded as a point mass m. Therefore, balloting of the projec- 
tile is not considered. A short time t after the projectile 
is fired, it lies at the barrel point x = £(t), as shown in 
Fig. 1. The axial velocity of the projectile relative to the 
barrel is v(t) = £ . The axial frictional force of the projectile 
on the barrel is FCt). The gas pressure driving the barrel is 
pCt). There is a tuning mass with constant coordinate n. 

Its center of mass is considered to lie on the axis of the 
barrel. Its moment of inertia about a transverse axis through 
its center of mass is 1^ . The breech is regarded as a rigid 
block with mass M and moment of inertia I about a transverse 
axis through its center of mass. The location of the center of 
mass is defined by two lengths, e and d shown in Fig. 1. The 
axial recoil displacement is u(t). The breech recoil spring 
and recoil dashpot are assumed to be linear. Their constants 
are Og , Pg) and their eccentricities are (e', e" ) , as shown 
in Fig . 1 . 

If there is initial bending of the barrel due to weight, 
unsymmetr ical thermal gradients, or manufacturing imperfections, 
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the bending is aggravated by inertial reaction between the barrel 
and the projectile. Axial friction F(t) also excites motion of 
the barrel. Another disturbance comes from the so-called 
"Bourdon" effect, which arises because, in the bent barrel, 
the area of the bore above the neutral plane, on which the gas 
pressure acts, is slightly greater than the area below the 
neutral plane. Still another effect contributing to the bending 
of the barrel is axial inertia of the recoiling barrel and breech. 

3. CONTRIBUTIONS OF VARIOUS PARTS OF THE SYSTEM AND OF VARIOUS 

EFFECTS TO THE EQUATIONS OF MOTION 

3.1. Strain Energy of the Barrel. The barrel may be 
regarded as a tapered elastic beam of length L. If bending 
strain energy only is included, the strain energy per unit length 
is proportional to the square of the curvature. If the effect 
of transverse shear deformation is included, an appropriate 
term for the shear energy must be added. In this report these 
energies are taken, respectively, as (1) 



U 



Bending 2 



U 



Shear 



El ( y" ) 2 dx 



GA(y ' -0 ) dx 



( 6 ) 



(7) 



where primes denote derivatives with respect to x, y = the 
lateral displacement of the barrel, 0 = the rotation of the 
barrel section (Note that if shear deformation is not included 
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the rotation 9 = y', the slope of the axis of the barrel), 

El = the bending rigidity and GA = the shear rigidity of the 

beam. By Hamilton's principle, these energies contribute 

terms to the components of generalized forces (Eq. 3) or 

equivalently to the stiffness matrix K (Eq. 4). 

In the usual method of finite elements (2), the barrel is 

divided into n intervals (elements) of length A A 0 , ..., A., 

l J 

..., A by points x-, , x 0 , ... x . , , where x is the 

n i ^ nti n + i 

coordinate of the muzzle (Fig. 1). These points need not be 
equally spaced. Then with y taken as a cubic polynomial in x 
in each interval of length A., and 9 taken as a corresponding 

Kj 

quadratic polynomial in x, the element stiffness matrix for the 
jth element is (see Appendix I), for j = 1, 2, ..., n, 



k. = 



E . I . 
J J 



A . (1 + r) . 
J J 



12 


6A 


-12 


6A 




6A 


A 2 (4+r ) 


-6A 


A 2 ( 2-r ) 


(8) 


-12 


-6 A 


12 


-6A 




6A 


A 2 ( 2-r ) 


-6A 


A 2 ( 4+r ) 


j 


\ . = 

J 


x j+l ~ : 


s . . A 
J 


ppropriate additions 



of the matrices k . are required to incorporate the element 

^ 3 

stiffnesses into the structural stiffness matrix K of Eq . (4) 

(2). The effects of shear deformation may be discarded by 
setting r = 0 in Eq. (8). 

3.2. Certain Effects Related to Pressure, Rigid Parts and 
Supports. The generalized coordinates are taken to be the dis- 
placements y., the rotations 9., the displacement £ of the 
J J 

projectile with respect to the breech, and the axial displacement 
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u of the barrel and the breech (as noted previously; see Fig. 1). 
The breech is regarded as a rigid block of mass M and moment of 
inertia T about its center of mass. The center of mass of the 
breech may be off center from the axis of the barrel. The 
eccentricity of the breech is specified by two lengths, d and e 
(Fig. 1). There is an axial spring with constant a_ and an 

D 

axial dashpot with constant at the breech. These elements 

are offset distances e f and e" from the axis of the barrel (Fig. 1) 

Movement of the breech transverse to the barrel is restrained 

by a spring with constant and a dashpot with constant y^. 

These are included among springs and dashpots with constants 

a. and y. that act at the nodal points on the barrel. Their 
3 3 

effects will be discussed later. 

3.2.1. Action of Gravity on the Breech. By Fig. 1, the 
descent of the center of mass of the breech due to displacement 
at x = 0 is 



u sina - y^ cosa + (d cosa - e sina)9^ 

where a is the angle of inclination of the barrel with respect 
to the horizontal. Therefore, the contribution of the action 
of gravity on the breech to the variation of the virtual work is 

Mg[(5u sina - Sy^osa + (d cosa - e sina)69 1 ] (9) 

By Eqs. ( 2 ) and (31, it is seen that the coefficients of 6u, 

Syf, <59^ in Eq. (_9) contribute to the right-hand side of Eqs. (3). 
With Eqs. (3) written in the form of Eqs. (4), the coefficients 
of <5u, 6 y , <59^ contribute to the right-hand side of Eqs. (4). 
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3.2.2. Contribution to <5W of Axial Force on Projectile. 

2 

The projectile is subjected to the axial force tra p-F, where p 



is the gas pressure behind the projectile, a is the radius of 
the bore, and F is the axial frictional force of the barrel on 
the projectile. The effect of the shock wave ahead of the 
projectile can be included in F. The contribution to 6W of the 
axial force on the projectile is, therefore, 



3.2.3. Contribution to 6W of Axial Movement of the Barrel. 
Since pressure p also acts on the breech, the contribution of 
the axial movement of the barrel to 6W is 



where M is the mass of the barrel, and where the term Mg(sina) 
is the contribution of the axial component of gravity. 

Accordingly, the net contribution to 6 W from the breech, 
the axial force acting on the projectile, the axial movement 
of the barrel, the breech spring, and the breech damper is, 
by Eqs . (9), (10), (11), 



The terms involving cig and ]jg in Eq. (15) come from the fact 

that the compression of the breech spring is u-e'9^. Consequently, 



(Tra^p-F ) ( 6 E, - 6u) 



( 10 ) 




( 11 ) 



SW = [-ag(u-e'9 1 )-yg(u-e M 9 1 ) 

+ (ii+M)g sinajSu +(ira^p-F)S£ 
-Mg(cosa)6y 1 + [Mg(d cosa - e sina) 
+a B e' (u-e ' 9^ )+ yge" (u-e"^ ) ] S0 1 



( 12 ) 
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the resisting forces of these elements are oig(u-e'0^) and 
y (u-e"0^). The coefficients of 6u, 6?, 6y^, 60^ contribute 
to the corresponding terms in Eqs . (3) or (4). 

3.2.4. Contribution to 6W of Dampers and Supporting 
Structure. The contribution to <5W from rotational and transla- 
tional dashpots at point x. is 

*J 

-y.y.6y. - yl 0.60. (13) 

3 3 3 3 J J 



where (y . , y '. ) are constants for the dashpots, and the dot 
3 J 

denotes time derivative. At a node where there is no dashpot , 




If the supporting structure can be regarded as a set of 
decoupled springs located at the nodes, the contribution of the 
springs at node j to 6W is 



-a .y . 6y . 
J J J 



a '. 8 . 6 0 . 
J J J 



(14) 



Accordingly, the contributions of the supports and dampers to 
6W is 



n+1 n+1 

j 

3=1 3=1 

Again, the coefficients of <5y . and 69. contribute in the appro- 

3 3 

priate manner to Eqs. (3) and (4). In particular, the coeffi- 
cients of y., 9. contribute to the matrix K and the coefficients 
3 3 — 

of y, 9. contribute to the matrix C in Eq. (4). 



£ 



( y '. 8 . +ct '. 0 . ) 6 9 . 
3 3 3 3 3 



(15) 



6W 



= - Z C Vj + V3 )6y 
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3.2.5. Action of Gravity on the Barrel. In the finite 
element method employed, the external loads (forces and couples) 
are applied at the nodes only. If a distributed load or a con- 
centrated load is applied between nodes, they must be replaced 
by equivalent systems of concentrated loads at the nodes. We 
require these equivalent loads to be derived in such a manner 
as to be consistent with the method used in deriving the stiff- 
ness matrix (Appendix A). Since the lateral component of the 
gravity load acting on an element of the barrel is uniformly 
distributed between nodes, the equivalent nodal force matrix 
for the jth element of the barrel is 



f 6 



f . 



B 



-m . A . 
_ J J 



g COSCt 
12 



1 

6 



-1 



j = 1,2, n (16) 



where m. = the mass of the jth element and A. = the length of 
3 5 

the jth element. By the method of finite elements, appropriate 

B 

additions of the matrices f . are required to incorporate the 

" w J 

element nodal forces into the driving force vector f(t), Eq . (4). 

The effect of the axial component of gravity on the barrel 
has been accounted for previously (see Eq. 11). 

3.2.6. Action of Gravity on the Tuning Mass. In general, 
the tuning mass is located at some point between nodes, Figs. 1 
and 2. Since the lateral component of gravity force on the 
tuning mass is m^ g cosa, this force must be transformed into 
equivalent nodal forces for the element in which the tuning mass 
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is located. By the method of finite elements (Appendix A), the 
equivalent nodal force matrix for the element in which the tuning 
element lies is 

(. 17 ) 

T 

where = the length of the tuning mass element and s^A^ = the 
distance from the left node of the element to the tuning mass. 
Equation (17) accounts for the lateral component of gravity 
acting on the tuning mass. In addition, the axial component of 
gravity contributes to SW in the form 

<5W = m T g (sina)Su (18) 

3.2.7. Action of Gravity on the Projectile. The absolute 
axial displacement of the projectile at time t is (*-u), Figs. 

1 and 3. Then, the contribution to oW of the axial component 
of gravity on the projectile is 

6W = mg (sina) C<5u - <$£) (19) 

The effect of the lateral component of gravity on the projectile 
is included in the same manner as for the tuning mass (Eq. 17) 
where A T , s T are replaced by A , s respectively. In addition, 
it is to be noted that s is a function of time t, since £ = £(t). 
In studying the effect of gravity on the projectile, it has been 
found by actual computation that this effect is extremely small 
for projectiles of the 20-40 mm class. 



- 3s 2 + 2s 3 



Cl 



- 2s + s )sA 









cosa 



3s 2 - 2s 3 



(_- s + s )sA 
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3.2.8. Axial Friction Effect of Projectile on Barrel. The 
axial frictional force F(t) exerted by the projectile on the 
barrel is assumed to be a known function of time t. Alternatively, 
it may be computed if the pressure acting on the projectile and 
the projectile acceleration are known as functions of time. 

Also, it may be expressed as F(£), where x = £(t) locates the 
projectile in the barrel. This is a distance s along the curve 
of the axis of the barrel (Fig. 4). 

There is a subtlety in the calculation of the work that 
force F performs on the barrel, because the projectile is 
sliding in the barrel. For example, if a grinding wheel acts 
on a fixed plate, the frictional force of the grinder performs 
no work on the plate because the plate does not move. However, 
the frictional force of the plate performs negative work on the 
grinder . 

The axial movement of the projectile during an infinitesimal 
time dt does not affect the work that the projectile performs 
on the barrel. Consequently, in computing 5W, we set dt = 0. 

The corresponding point on the curve y + 6y consequently also 
lies a distance s along the curve. The component of F along the 
virtual displacement determines 6W (Fig. 4). We have 




/I + (y ' ) 2 dx 



o 



Also, along the curve y + 6y, 




o 



22 



Since y < < 1 , 






( y ' ) 2 Cbc 



and 



x r 



x 2 + 2 



(y ’ 2 + 2y ' 6y ’ ) dx 



Therefore, 



X - 



x r 



X 2 X 1 






( y ’ ) 2 dx 



i _ 

2 



( y ’ 2 + 2y ’Sy’ ) dx 



and since x^ = £(t), 



C +£ 



£ = y’ dx - ^ / ( y ' + 2y ' <5y ' ) dx 



and we have 



1 / „ , 2 



,2 



e = o y' dx - ^ / ( y ’ + 2y ' 5y ' ) dx 



C+e 



/ 



, 2 , 
y dx 
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or 



?+£ 



i 



£ = - / y ' 6y 'dx - ^ 



/ 



y' 2 dx 



However, since e is small, 



C+£ 



/ 



y’ 2 dx = e [y ' ( ? ) ] 2 



Hence , 



C 

S 



£ { i + J[y'(C)] 2 } = -/ y’Sy'dx 



-I 



and since y' << 1 



e » -J y ' 6 y ' dx 
o 



Using the trapezoid rule to evaluate this integral, and setting 
C = x r , we get 

£ = - ^-[ ( 0 -^ 5 + 02 + ( + ®3 

+ ( ^4 )^3 + • • • • 

+ (0 ,60 , +0 60 )A , ] 

v r-1 r-1 r r r-1 

The contribution to 5W from axial friction is the component 
of F in the x direction times £ plus the component of F in the 
y direction times 5y . Thus, we have approximately 

6W = Fe + F0 ay 
r J r 
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